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EXECUTIVE  SUMMARY 


APPLICATION  OF  SENSITIVITY  VECTORS  TO  THE  MEASUREMENT 
AND  MODELLING  OF  MAGNETOSTATIC  FIELDS 

OBJECTIVE 

The  objective  of  this  report  is  to  develop  the  analytic  techniques 
required  for  quantitative  assessment  of  the  sensitivity  of  magnetic  field 
measurements  to  quasi-static  magnetic  sources. 

APPROACH 

To  accomplish  this  objective  we  first  use  the  concept  of  the  gradient 
vector  to  define  a  "sensitivity  vector"  for  a  given  field  measurement,  and 
then  relate  this  vector  to  the  results  obtained  from  the  reciprocity  theorem 
of  electromagnetism.  The  derivations  will  then  be  extended  from  magnetic 
dipoles  to  higher  order  moments.  Several  non-linear  cases  will  be  examined, 
and  techniques  will  be  developed  to  account  for  the  presence  of  noise. 
RESULTS 

In  this  report,  the  mathematical  basis  of  sensitivity  vectors  has  been 
developed,  with  specific  application  to  the  measurement  and  modelling  of 
magnetic  fields.  In  the  course  of  analyzing  several  examples,  the  following 
questions  were  answered: 

1.  Given  two  magnetic  field  measurements,  what  criteria  must  be 
satisfied  to  allow  determination  of  the  two  dipole  components 
consistent  with  that  field? 

2.  For  a  measurement  of  the  vector  magnetic  field  at  a  fixed 
distance  from  a  dipole  source,  what  is  the  optimum  magnet¬ 
ometer  position? 

3.  How  do  various  configurations  of  single  axis  magnetometers 
compare  in  terms  of  their  ability  to  determine  model  parameters, 
and  how  can  this  be  assessed  quantitatively? 

4.  How  can  the  presence  of  noise  be  included  in  such  an  analysis? 


-  vi  - 


5.  What  is  the  relationship  of  the  sensitivity  vector  concept  to 
the  reciprocity  theorem  of  electromagnetic  fields? 

6.  How  can  sensitivity  vectors  be  used  to  study  multipole  models? 

7.  What  is  the  interpretation  of  a  sensitivity  vector  for  a 
non-linear  model? 

8.  What  is  the  position  dependence  of  a  vector  magnetometer, 
and  how  does  it  vary  with  position? 

9.  Where  can  a  magnetometer  be  placed  relative  to  the  model  to 
insure  adequate  signal  to  noise  ratio? 

Based  on  the  ability  of  the  sensitivity  vector  approach  to  answer  these 

questions  quantitatively,  it  appears  that  this  type  of  analysis  may  be 

valuable  for  optimizing  magnetometer  array  configurations.  It  is  reassuring 

that  the  results  obtained  in  this  report  are  consistent  with  both  a  modeller's 

intuition  and  more  abstract  mathematical  analysis,  particularly  in  that  this 

method  can  be  readily  extended  to  more  complicated  systems  where  intuition 

fails. 

RECOMMENDATIONS 

Two  specific  recommendations  follow  from  these  results: 

1)  Interactive  computer  code  should  be  developed  to  allow  accurate 
and  rapid  analysis  of  magnetometer  sensitivity  and  comparison  of 
magnetometer  configurations. 

2)  The  sensitivity  vector  analysis  should  be  extended  to  include  the 
magnetic  field  from  electric  current  distributions,  and  possibly 
the  electric  field  from  these  currents. 
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I .  INTRODUCTION 


The  purpose  of  this  report  Is  to  develop  the  analytic  techniques  re¬ 
quired  for  quantitative  assessment  of  the  sensitivity  of  magnetic  field 
measurements  to  quasi-static  magnetic  sources.  To  do  this,  we  will  first 
use  the  concept  of  the  gradient  vector  to  define  a  "sensitivity  vector" 
for  a  given  field  measurement,  and  then  relate  this  vector  to  the  results 
obtained  from  the  reciprocity  theorem  of  electromagnetism.  The  deriva¬ 
tions  will  then  be  extended  from  magnetic  dipoles  to  higher  order  moments. 
Several  non-linear  cases  will  be  examined,  and  techniques  will  be  deve¬ 
loped  to  account  for  the  presence  of  noise. 
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II.  THE  CONCEPT  OF  A  SENSITIVITY  VECTOR 
Let  us  consider  an  object  that  is  a  source  of  electric  and  magnetic 
fields.  We  are  able  to  make  electromagnetic  measurements  at  some  distance 
from  the  object,  and  wish  to  use  these  measurements  to  obtain  a  mathema¬ 
tical  description  of  the  field  sources.  We  can  accomplish  this  by 
defining  a  hypothetical  model  for  the  sources  and  by  adjusting  various 
model  parameters  until  the  fields  produced  by  the  model  match  the  observed 
ones.  In  practice,  this  process  is  complicated  by  linear  dependence  of 
various  model  parameters  for  a  particular  set  of  measurements,  by  the 
presence  of  noise  in  the  measured  data,  and  by  the  inability  of  the  model 
to  explain  certain  details  of  the  fields.  The  last  one  of  these  compli¬ 
cations  can  be  remedied  only  by  altering  the  model  and  will  not  be  consi¬ 
dered  in  this  report.  The  first  two,  linear  dependence  and  noise,  can  be 
addressed  using  sensitivity  vectors. 

* 

We  will  treat  each  measurement  as  a  scalar  .  Measurement  of  a  vector 
magnetic  field  is  equivalent  to  three  single  axis  field  measurements  that 
determine  three  orthogonal  vector  components.  Other  scalars  that  we 
might  measure  are  the  component  of  the  magnetic  field  parallel  to  the 
earth'p  geomagnetic  field,  B  Bearth  ,  the  field  magnitude  [b|,  or  electric 

field  components  E^  and  magnitude  | E | .  In  general,  we  will  make  field 
measurements  at  points  r^.  An  arbitrary  number  of  different  measure¬ 
ments  can  be  made  at  a  single  point,  in  which  case  F  =  F  but  r  =  r 

*  J  1  J  • 

The  process  of  modelling  involves  using  these  measurements  to  specify  the 
parameters  of  a  model  which  reproduces  the  fields  to  the  desired  accuracy. 
Suppose  the  model  is  located  at  r1  and  is  specified  by  n  model  parameters 
Mj .  As  an  example,  a  magnetic  dipole  model  at  a  known  location  has  three 
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*i.e.,  a  scalar  quantity  as  opposed  to  a  vector  quantity,  not  to  be  con 
fused  with  a  total  field  magnetometer. 
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parameters  =  m^,  M2  =  m^  and  =  in^.  Addition  of  a  quadrupole  adds 
five  more  parameters  to  the  model. 

If  111  measurements  are  made  to  determine  the  n  parameters  M^ ,  we 
can  say  that  the  field  or  measurement  space  has  n  dimensions  while  the 
source  or  model  space  has  m  dimensions.  If  the  F's  and  M's  are  linearly 
related,  we  can  write 


T1 . T1M  \  /  M1 


TN1  *  *  •  "  ^NM  /  V’m 


(2.1) 


which  in  matrix  notation  is 


*  - 


(2.2) 


where  f  is  a  n  x  1  matrix;  ill  is  a  m  x  1  matrix,  and  ^  is  the  m  x  n  trans¬ 
fer  matrix.  If  the  F's  and  M's  are  not  linearly  related,  as  would  be  the 
case  if  dipole  location  was  a  parameter  in  the  model,  we  can  linearize 
the  equations  about  a  point  in  model  space.  Non-linear  examples  will  be 
treated  in  a  later  section. 

In  particular,  each  element  T  represents  the  value  of  the  ith 
field  measurement,  i.e.  F^,  if  only  one  model  term  M  contributes,  i.e. 

‘*k  “  ^kj*  lhus  a  horizontal  row  of  ^  corresponds  to  the  field  equation 
for  tiie  complete  model  and  describes  how  a  particular  measurement  is 
affected  by  each  source  term,  while  each  column  of  ^  describes  how  a  par¬ 
ticular  source  term  affects  each  of  the  measurements.  If  there  are  more 
measurements  than  model  parameters,  Eq.  (2.1)  will  represent  an  over¬ 
determined  set  of  equations.  However,  it  is  possible  that  a  pair  of  model 
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parameters  may  not  be  totally  independent,  resulting  in  coupled  columns 
in  T.  Then  T  might  be  singular  so  that  the  set  of  equations  will 

*\j  *\i 

not  have  a  unique  solution  for  M  given  T  and  F.  Similarly,  if  two  or 

% 

more  field  measurements  are  linearly  related,  two  or  more  rows  of  T  would 
be  linearly  dependent.  This  will  pose  a  problem  if  there  is  an  insuffi¬ 
cient  number  of  independent  field  measurements  to  solve  Eq.  (2.1)  for  the 
model  parameters. 

As  a  simple,  two  dimensional  example  of  how  a  set  of  data  relates  to 
a  model,  let  the  model  be  a  magnetic  dipole  m  located  at  a  known  point  r  . 
We  will  attempt  to  make  measurements  of  at  points  r^  and  r ^  in  order 
to  determine  the  unknown  dipole  components  m^  and  m^.  The  magnetic 

induction  B  is  given  by 

'  .  V ' 


(2.3) 


so  that  B^  satisfies  the  equation 


B  (r)  =  7^ 
x  4  it 


[3  -  »•)  +  ■„<?  -  y'-]  (x  _  _[f) 


r  -  r 


•  !  3 


(2.4) 


For  two  B^  measurements,  we  will  have  two  simultaneous  equations,  linear 
in  m,  that  are  of  the  form  of  Eq.  (2.4).  We  must  determine  the  values  of 
and  for  which  these  two  equations  have  a  unique  solution.  As  we 

will  see,  sensitivity  vectors  will  allow  us  to  do  this  and  also  to  identify 
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the  optimum  values  of  and  Tj. 

We  can  define  vector  spaces  for  both  the  field  measurements  and  the 
model  parameters,  so  that  the  elements  of  ^  and  ?(  become  the  components 
of  vectors  ?  and  il.  As  an  aid  in  understanding  the  relationship  of 
several  measurements  to  the  model,  we  can  define  a  gradient  operator  in 
model  space  by 


^’3^i+3^j  +  3M^k+‘ 


*  *3M 


(2.5) 


where  i,  j,  k,  m  are  the  unit  vectors  in  each  direction  in  model  space. 
In  our  two-dimensional  example, 


7  -  -i-  i  +  -L-  j 

M  3m  3m  ■* 

x  y 


(2.6) 


The  sensitivity  vector  for  measurement  F^  is  defined  by 


5i  ■  5H  fl 


(2.7) 


The  interpretation  of  is  straightforward:  indicates  the  direction 

of  the  change  in  model  ?(  that  produces  the  maximum  change  in  the  measure¬ 
ment  F^.  Any  change  in  ?f  that  is  perpendicular  to  will  be  undetected 
by  F^.  Thus,  the  direction  of  the  sensitivity  vector  can  be  used  to 
determine  the  source  configuration  for  which  has  the  maximum  sensiti¬ 
vity  and  the  configuration  which  is  not  detected  by  F  .  The  magnitude  of 
the  sensitivity  vector  is  proportional  to  the  maximum  field  produced  at 
the  measurement  point  by  a  given  strength  source  and  can  be  used  to 
determine  the  relative  sensitivities  of  several  measurements. 

For  linear  systems,  the  components  of  the  sensitivity  vector  for  a 
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particular  measurement  is  in  fact  a  row  of  T.  To  show  this,  we  can 


write  Eq.  (2.2)  as 


Fi  =  j  TijMj 


(2.8) 


and  Eq,  (2.7)  as 


m  m  j  n  » 

S  <Fi>Sk'  S  “k(?iTy  Mj) 


(2.9) 


where  ei,  is  the  kth  unit  vector. 


Since  the  T. .  are  independent  of  M,  this  becomes 


(2.10) 


j,k=l  K 


^  Tij  6jkek  ^  TikSk 

j>k=l  J  J  k=l 


"►  th  ^ 

Thus  is  a  vector  in  model  space  whose  components  are  the  i  row  of  T. 

While  we  have  obtained  this  simple  result  from  an  apparently  circuitous 

route,  it  provides  us  with  a  better  understanding  of  the  significance  of 

each  row:  and  the  iC^  of  ^  can  be  used  to  determine  the  sensitivity  of 

to  changes  in  the  model  M,  and  vectors  can  be  plotted  to  help  visualize  how 
measurement  sensitivity  depends  upon  measurement  location.  For  non-linear 
systems,  which  can  not  be  described  using  equations  of  the  form  of  Eq-.  (2.1) 
and  (2.2),  we  will  show  subsequently  that  the  sensitivity  vectors  are  still 
valid  and  can  be  readily  determined  analytically  or  numerically  without 
explicitly  linearizing  the  system  of  equations. 


III.  A  TWO-DIMENSIONAL  DIPOLE  MODEL 


In  the  simple  two-dimensional  model  of  a  magnetic  dipole  m  with  com¬ 
ponents  m  and  m  ,  the  field  in  the  x-direction  is  given  by  Eq.  (2.4). 
x  y 

The  sensitivity  vector  becomes 


■>  v  ■>  -f 

S(r ,  r')  =  V  B  (r) 

m  X 


(3.1) 


or,  after  dropping  the  ^  for  convenience, 

S(r ,  ?’)  =  —  ■■ -  {[3(x  -  x')2  - 

I  r  -  r*  I  5  1 


r  -  r 


I2]i 


+  3(x  -  x')(y 


-  y').i} 


(3.2) 


If  the  dipole  is  at  the  origin,  r'  =  0  and 


S(r,  r)  =  — —  |(3x2  -  r2)i  +  3xyj| 
r5  " 


(3.3) 


We  can  use  Eq.  (3.3)  to  calculate  the  sensitivity  vectors  for  measurements 


of  Bx  in  the  xy  plane.  It  follows  that 


—  (3x2  -  r2) 

T  ^ 


=  3xy 


(3.4) 


If  ali  measurements  are  made  at  the  same  distance  from  origin  then 

1 


|r|  «  ^ x2  +  y2  =1  and 


S  =  3x2  -  1 


S  =  3xy 

y 


(3.5) 


The  sensitivity  vector  components  and  magnitudes  at  sixteen  points  are 

listed  in  Table  3.1  and  plotted  in  Figure  3.1.  Several  important  features 

are  immediately  obvious.  At  points  A  and  E,  a  measurement  of  B  detects 

x 
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fields  from  m  but  not  from  m  .  The  sensitivity  at  A  is  twice  that  at  E. 

x  y 

At  point  C,  niy  is  detected  but  not  n>x.  The  sensitivity  components  Sx  and 
Sy  are  equal  to  each  other  at  point  B  and  at  point  0,  but  the  magnitude 

of  3  is  1.81  at  B  and  only  1.10  at  D. 

If  Bx  is  to  be  measured  at  two  locations  for  the  purpose  of  deter¬ 
mining  both  mx  and  m^,  the  figure  can  be  used  to  determine  which  combina¬ 
tions  will  not  provide  a  solution  for  both  m  and  m  .  If  the  two  B 

r  x  y  x 

measurements  are  made  at  the  two  A  points,  or  at  the  two  E  points,  it 


will  not  be  possible  to  determine  m^  because  both  sensitivity  vectors  are 
parallel  to  the  mx  axis.  Similarly  a  measurement  at  A  and  another  at  E 


will  not  work  Any  two  C  measurements  have  parallel  sensitivity  vectors 
and  are  thus  unsatisfactory.  These  are  measurement  combinations  for 
which  the  matrix  'i1  in  Eq.  (2.2)  has  no  inverse. 

The  figure  can  also  be  used  to  determine  which  combinations  provide 
a  good  determination  of  n»x  and  m^.  Two  parallel  vectors  were  shown  to  be 
bad.  Two  perpendicular  vectors  may  be  good,  but  will  two  larger,  non¬ 
perpendicular  vectors  be  better?  A  measure  of  the  combined  "information" 
contained  in  two  sensitivity  vectors  can  be  defined  in  terms  of  the 
magnitude  of  their  cross-product 

v  “  \i±  X  $2 I  (3.6) 

If  ^  and  ^  are  parallel,  V  ■  0.  If  and  ^  are  perpendicular,  V 


I 

0 

I 

1! 

3 

.1 

il 


l 


In  the  two-dimensional  example 
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For  two  B  measurements  in  adjacent  quadrants,  V  »  3.28,  for  two  D  measure¬ 
ments,  V  -  1.22,  implying  that  in  the  presence  of  noise,  it  would  be 
significantly  better  to  measure  at  adjacent  B  points  than  at  adjacent 
D  points.  The  values  of  V  for  pairs  of  measurement  points  in  the  first 
and  second  quadrants  are  listed  in  Table  3.2.  As  might  be  expected  from 
the  figure,  adjacent  D  and  E  points  have  a  low  V  value  of  0.78.  Note  the 
symmetry  of  the  table  across  the  heavy-lined  boxes. 

As  a  second  example,  suppose  we  have  a  single  vector  magnetometer 

that  measures  B  and  B  at  a  single  point  (x,  y) .  If  the  model  is  a  dipole 
x  y 

located  at  the  origin  and  lying  in  the  x-y  plane,  the  sensitivity  vector 
for  each  field  component  is  given  by 


S(Bx)  =  ~  £(3x2  -  r2)i  +  (3xy)j] 

S(By)  -  [(3xy)i  +  (3y2  -  r2)j] 


(3.8) 


The  cross  product  V  becomes 


v  -  IVBx>VV‘  Sj,(Bx)Sx(By)| 

-  |<3x2  -  r2)(3v2  -  r2)  -  (3xy)(3xy)|  *  (3.9) 

r  r 

This  result  indicates  that  the  ability  of  a  vector  magnetometer  to  deter¬ 
mine  the  components  of  the  dipole  depends  solely  on  the  distance  of  the 
magnetometer  from  the  dipole,  and  not  upon  the  relative  orientation  of 
the  dipole  and  the  magnetometer.  This  is  apparent  in  Fig.  3.2,  which 


FIRST  QUADRANT 


_..,i  I  II*  w  -nfT  riT  -  --  Hipmwpi'  Win  mush 


Figure  3.2 .  The  sensitivity  vectors  for  measurements  of  B  and  B  at  unit 


distance  from  a  fixed  magnetic  dipole.  The  coordinates  of  the  points  are 


listed  in  Table  3. 2, 
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shows  the  Bx  and  B^  sensitivity  vectors  for  a  vector  magnetometer  at  a 
constant  distance  from  the  dipole.  At  locations  A  and  E,  the  B^  magneto¬ 
meter  is  sensitive  only  to  m  ,  while  the  B  magnetometer  is  sensitive 

x  y 

only  to  ra  .  Thus  the  calculation  of  m  from  B  can  be  performed  with  equal 

y 

accuracy  at  all  points  on  the  circle  in  the  figure.  We  will  show  later 
that  the  location  of  a  vector  magnetometer  becomes  significant  when  using 
a  model  with  higher-order  moments. 
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IV.  THE  INFORMATION  MATRIX 

While  the  quantity  V  provides  a  measure  of  sensitivity  and  the 
ability  to  invert  the  transfer  matrix  to  determine  the  model  parameters, 
we  need  more  complex  models  and  thus  must  find  a  multidimensional  equi¬ 
valent  of  V.  Extension  of  Eq.  3.5  to  three  dimensions  is  straightforward, 
with  V  defined  by  the  vector  triple  product 


V  -  •  (s2  x  s3) 


S2  •  (S3  x  -  S3  •  (S3  x  S2)  (4.1) 


Note  that  V  is  equal  to  the  volume  in  model  space  enclosed  by  the  three  sensitivity 
vectors.  If  any  two  of  the  vectors  are  parallel,  V  *  0.  We  see  that  for 
a  three  parameter  model,  V  measures  how  well  the  sensitivity  vectors  for 
a  set  of  measurements  "span"  model  space.  In  the  absence  of  noise  or 
computational  inaccuracy,  an  n  parameter  model  can  ideally  be  specified 
using  n  measurements  that  have  orthogonal  sensitivity  vectors  spanning 
model  space.  In  the  three-dimensional  case,  three  measurements  are  re¬ 
quired  but  additional  measurements  may  increase  V  either  because  some  of 
the  measurements  may  have  low  sensitivity  or  non-orthogonal  sensitivity 
vectors. 

In  order  to  extend  V  to  more  than  three  dimensions,  we  need  to  intro¬ 
duce  the  information  matrix.  The  component  of  the  sensitivity  vector 
for  the  ith  measurement  is  (3^  and  is  given  by  Eq.  (2.10) 

"  Tij  <4'2> 

so  we  can  rewrite  Eq.  (2.8)  as 


!i  ■  j  Tij  Mj  ■  «*t>j  Hj  ■  V » 


(4.3) 
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A.  T  “► 

where  is  the  transpose  of  the  column  matrix  equivalent  to  S^.  We 

can  multiply  Eq.  (4.3)  on  the  left  by  to  obtain,  for  an  m  dimensional 

model , 


=  *  * T  a 


(4.4) 


where  is  termed  the  "information  matrix"  and  is  an  m  x  m  square,  symmetric 
matrix  particular  to  the  ith  measurement.  If  there  are  n  measurements, 
the  n  equations  of  this  form  can  be  summed  to  yield 


£  <W>a  ■ 


i=l 


n 


a 


(4.5) 


which  can  be  solved  for  ft 


ft 


i-l 


-1  n 


E*J  Ev^Ev* 


i-l 


i=l 


(4.6) 


with  K  being  the  m  x  m  information  matrix  for  the  entire  set  of  measure¬ 
ments.  Equation  (4.6)  represents  a  solution  to  the  "inverse  problem"  in 
which  field  data  are  used  to  determine  certain  source  parameters.  In  the 
equation,  the  n  F  will  be  known,  since  these  are  the  field  measurements, 
and  the  and  R  can  be  calculated  for  the  chosen  model.  However,  if  ft 
does  not  have  an  inverse,  we  cannot  determine  ft. 

As  an  example  of  this,  let  us  return  to  the  two-dimensional  dipole 
example  in  Section  3.  If  measurements  are  made  at  points  A  and  E,  the 
sensitivity  vectors  and  information  matrices  will  be 


i 
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^  1 
R1  “  S1S1 


0.0 


0.0 


(4.7) 


% 

R 


0.0 


0.0 


Because  R  has  a  zero  determinant,  it  has  no  inverse  and  we  can  not  deter¬ 
mine  both  dipole  components.  This  is  as  expected,  since  for  these  two 
measurements,  V  ■  0. 

We  can  make  the  connection  between  the  determinant  of  ft  and  V  by 
writing 


and  computing  det 


1*1 


det |ft| 


<S21S12  -  S11S22>2  '  V2 


(4.8) 


(4.9) 


to  find  that  our  n-dimensional  V  is  simply 

V  -  (detjftl)*4 


(4.10) 


We  now  see  that  V,  which  provides  a  measure  of  how  well  a  set  of  measure¬ 
ments  span  the  model  space,  is  also  a  measure  of  how  readily  the  informa¬ 
tion  matrix  can  be  inverted  as  required  to  solve  the  inverse  problem.  Given 
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%  'V 

this  introduction  to  S^,  V,  and  R,  we  are  now  prepared  to  work  several 
examples. 

4.1.  The  Two-Dimensional  Dipole 

Equation  4.4  ana  the  data  listed  in  Table  3.1  can  be  used  to  compute 
the  information  matrix  for  each  measurement  in  quadrants  I  and  II  in  Fig. 

3.1,  and  the  information  matrix  for  all  pairs  of  measurements. 

The  results,  listed  in  Table  4.1,  are  consistent  with  Table  3.2  and  show 
that  certain  pairs  such  as  A1-A2,  Al-El  and  C1-C2  have  a  singular  infor¬ 
mation  matrix.  The  pair  B1-B2  had  the  largest  value  V  in  Table  3.2  and 
is  now  seen  to  have  a  diagonal  information  matrix  with  large,  equal 
eigenvalues.  In  this  case,  R  can  be  written  as 


2f  -  A  * 

where  A  is  the  eigenvalue  and  1  is  the  identity  matrix, 
reduces  to 

A  -  p 

i=l  1 

th  'V 

The  j  '  component  of  M  becomes 


(4.11) 

Equation  4.6 

(4.12) 


■  1-1  5‘Vj  i  ■  £  t4J 


(4.13) 


We  recognize  the  summation  as  a  matrix  multiplication  and  write 


M  =  A"1  *T  £ 


(4.14) 


•hus  for  properly  chosen  measurement  points,  solution  of  the  inverse 
problem  reduces  to  a  trivial  matrix  multiplication  using  the  transpose  of 
tne  original  ^  matrix.  This  is  the  motivation  for  finding  sets  of 
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measurement  points  where  the  sensitivity  vectors  are  orthonormal. 

As  an  example  of  this,  let  us  assume  that  we  can  make  one  B^  field 

measurement  at  the  arbitrary  point  (0.75,  0.66),  between  points  B  and  C  in 

a. 

Fig.  3.1,  and  that  we  want  to  find  the  second  point  so  that  R  is  diagonal 

%  % 

and  has  equal  eigenvalues.  Using  Eq.  (3.3),  we  find  that  and  will 
be 


-V  /o.69\ 


The  information  matrix  for  these  two  measurements  is 


(4.15) 


<v 

R  = 


(2x2  -  y2)2  +  0.48r10 
3xy(2x2  -  y2)  +  1.03r10 


3xy(2x2  -  y2)  +  1.03r 
(3xy)2  +  2. 22r 1  0 


(4.16) 


Rather  than  use  Eq.  (4.16)  to  find  the  values  of  x  and  y  where  R  is  di¬ 
agonal  and  has  equal  eigenvalues,  it  is  more  convenient  to  determine  the 
points  where 


Itj  *  |S2I  and  Si  •  S2  =  0 


(4.16) 


From  Eqs.  (4J.5) ,  it  follows  that  Eq.  (4.16)  is  satisfied  by 


52  -  +1.49  x  +0.69  y 


(4.17) 


Visual  examination  of  Fig.  3.1  shows  that  the  points  with  this  sensitivity 
will  be  located  approximately  as  follows 

Quadrants  I  &  Ill  between  0  &  E,  r  <  1 
Quadrants  II  &  IV  between  A  &  B,  r  >  1 


(4.18) 


- 

' 
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We  can  now  proceed  to  calculate  the  locations  exactly  using  Eqs.  (4.15) 
and  (4.17) 


(2x2  -  y2)  =  +  1.49 


(4.19) 


(3xy) 


0.69 


The  most  direct  method  of  solving  these  simultaneous,  non-linear  equations 
is  to  solve  first  for  r  and  then  for  x  and  y,  which  gives  the  values 
r  -  1.057  r  =  0.856 

x  =  t  1.013  x  -  ±  0.125  (4.20) 

y  -  +  0.299  y  «  t  0.847 

These  data  are  plotted  in  Fig.  4.1.  The  significance  of  this  calculation 
is  that  the  results  were  predicted  approximately  using  Fig.  3.1.  Thus  one 
strength  of  sensitivity  vectors  is  their  suitability  visualizing  measure¬ 
ment  sensitivity. 


4.2  The  Single  Vector  Magnetometer 

We  can  now  use  the  information  matrix  to  examine  the  vector  magnetometer 

example  discussed  in  Section  3.  If  the  magnetometer  measures  B  and  B  at  a 

x  y 

point  (x,  y) ,  Eqs.  (3.8)  can  be  used  to  compute  the  information  matrices  for 

measurements  of  B  and  B  : 

x  y 

%  'V  vp 

Ri  "  si  si 


-v  ,  I  (3x2  -  r2)2 

R  -  — 

x  r10  V  3xy(3x2  -  r2) 


(4.21) 


**<■*•'  t 
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((3xy)2  3xy(3y2  -  r2) 

3xy(3y2  -  r2)  (3y2  -  r2)2 


(4.21  cont'd) 


Since  V  is  determined  by  r  but  is  independent  of  x  and  y,  this  shows  that 
there  is  no  optimum  position  for  the  vector  magnetometer  for  a  fixed  dis- 
tance  from  the  dipole.  Equation  (4.22)  also  shows  that  R  is  diagonal  for 
measurements  on  either  the  x  or  y  axes. 

In  a  later  section,  we  will  extend  this  example  to  include  different 
noise  in  the  x  and  y  axis  magnetometers. 

4.3  Multiple  B  and  B  Magnetometers 

_ 5 _ JL _ 

The  analysis  techniques  developed  in  the  preceeding  sections  can  be 
applied  to  the  row  of  either  11  B^  or  11  single-axis  magnetometers  shown 
in  Fig.  4.2.  The  sensitivity  vectors  and  information  matrices  for  each  mag¬ 
netometer  are  listed  in  Tables  4.2  and  4.3,  with  the  source  assumed  to  be 
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TABLE  4.2 

Sensitivity  vectors  and  information  matrices  for  the  line  of  magneto¬ 
meters  in  Figure  4.2.  S  ■  (3x2  -  r2)r-5,  S  ■  3xy  r-5 

x  y 

r  S  q  ■>  ^ 

Position  Coordinates  x  y  |s|  R 


-4.0,  1.0  4.12  0.026  -0.010  0.028 


-3.0,  1.0  3.16  0.054  -0.028  0.061 


-2.0,  1.0  2.24  0.125  -0.107  0.165 


-1.0,  1.0  1.41  0.177  -0.530  0.559 


-0.5,  1  1.12  -0.286  -0.859  0.905 


0.0,  1.0  1.0  -1.000  0 


1.000 


0.5,  1.0  1.12  -0.286  0.859  0.905 


1.0,  1.0  1.41  0.177  0.530  0.559 


2.0,  1.0  2.24  0.125  0.107  0.165 


3.0,  1.0  3.16  0.054  0.028  0.061 


4.0,  1.0  4.12  0.026  0.010  0.028 


0.001 
0.000 
0.003 
^-0. 002 
0.016 
v -0.013 
0.031 
v-0.094 
0.082 
0.246 
1.000 
0.000 
'  0.082 
^-0.246 
0.031 
0.094 
0.016 
0.013 
0.003 
0.002 
0.001 
0.000 


0.000 

0.000 

-0.002 

0.001 

-0.013 

0.011 

-0.094 

0.281 

0.246 

0.738 

0.000 

0.000 

-0.246 

0.738 

0.094 

0.281 

0.013 

0.011 

0.002 

0.001 

0.000 

0.000 


TABLE  4.3 

Sensitivity  vectors  and  information  matrices  for  the  line  of  magnetomers 

in  Figure  4.2.  S  =  3xy  r“5,  S  -  (3y2  -r2)r-5. 
x  y 


Position  Coordinates 


r  S 


-4.0,  1.0 


-0.010  -0.012  0.015 


0.000  0.000 


0.000  0.000 


-3.0,  1.0 


-0.028  -0.022  0.036 


0.001  0.001 


0.001  0.000 


-2.0,  1.0 


-0.107  -0.036  0.113 


0.011  0.004 


0.004  0.001 


-1.0,  1.0 


-0.530  0.177  0.559 


0.281  -0.094 


•0.094  0.031 


-0.5,  1.0 


-0.859  1.002  1.319 


0.738  -0.861 


-0.861  1.004 


0.0,  1.0 


2.000  2.000 


0.000  0.000 


0.000  4.000 


0.5,  1.0 


0.859  1.002  1.319 


0.738  0.861 


0.861  1.004 


1.0,  1.0 


0.530  0.177  0.559 


0.281  0.094 


0.094  0.031 


2.0,  1.0 


0.107  -0.036  0.113 


0.011  -0.004 


,-0.004  0.001 


3.0,  1.0 


0.028  -0.022  0.036 


0.001  -0.001 

-0.001  0.000 


4.0,  1.0 


0.010  -0.012  0.015 


0.000  0.000 


0.000  0.000 


I 
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a  dipole  located  at  the  origin.  The  sensitivity  vectors  are  plotted  in 
Figs.  4.3  and  4.4  both  as  scalar  functions  of  x  and  as  vectors  at  each 
measurement  point. 

The  plots  of  S  and  S  versus  x  in  Fig.  4.3  can  be  seen  to  correspond 
r  x  y 

I  , 

to  plots  of  the  field  of  the  mx  and  my  dipole  components,  respectively. 

This  will  be  explained  in  terms  of  the  reciprocity  theorem  in  Section  6.  i 

Examination  of  the  S  and  S  curves  can  be  used  to  choose  optimum  loca- 

x  y  i  i 

tions  for  a  Bx  magnetometer.  For  example  Bx  measured  at  F.  and  G  would  be  / 

good  for  determining  m  ,  but  Bx  at  F  would  be  more  sensitive  to  mx»  ,  I  I 


Choice  of  measurement  locations  is  simplified  by  the  vector  plot  in  the 


and  appear  to  span  the  dipole  space  relatively  well.  Addition  of  points 


d  and  H  will  make  only  small  improvements,  since  these  vectors  are  sign! 


ficantly  smaller  than  the  parallel  vectors  at  G  and  E,  respectively 


The  plots  in  Fig.  4.4  show  that  the  sensitivity  of  B  to  m  is  the 


Thus,  for  a  row  of  magnetometers  displaced  from  the  source  in  y 


direction,  measurements  of  B  are  superior  to  measurements  of  B  .  The 


vector  plot  in  Fig.  4.4  shows  that  B  measurements  at  locations  E  and  G 


have  large,  nearly  perpendicular  sensitivity  vectors.  Both  figures  show 


that  there  is  a  serious  loss  in  sensitivity  for  measurement  locations  for 


which  x  >  2y 


Given  the  sensitivity  vector  data  and  information  matrices  listed  in 


Tables  4.2  an(j  4.3,  we  can  calculate  the  quantity  V  for  various  combina' 


tions  of  B  and  B  measurements,  as  listed  in  Table  4.4.  Several 


important  points  become  apparent  in  the  table.  A  vector  magnetometer  at 


(0.0,  1.0)  has  a  V  of  2.0,  greater  than  any  of  the  combinations  of  2  to  11 


psw*«*kw** 
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TABLE  4.4 

The  quantity  V  for  various  combinations  of  B  and  B  magnetometers.  The 

x  y 

locations  are  shown  in  Fig.  4.2.  The  mean  distance  of  the  measurements 
to  the  dipole  is  listed  under  r. 


Locations 

r 

B 

X 

B 

y 

B  &  B 
x  y 

Plot  1 

A 

4.12 

0 

0 

0 

B 

3.16 

0 

0 

0.002 

C 

2.24 

0 

0 

0.016 

0 

1.41 

0 

0 

0.249 

2 

E 

1.12 

0 

0 

1.025 

3 

F 

1.0 

0 

0 

2.000 

4 

E,  G 

1.12 

0.492 

1.722 

2.390 

5 

E,  F,  G 

1.08 

1.311 

2.978 

4.445 

6 

D,  F,  H 

1.27 

0.773 

1.511 

2.740 

7 

D,  E,  F, 

G, 

H 

1.21 

1.581 

3.517 

5.144 

8 

F,  G 

1.06 

0.859 

1.718 

3.174 

9 

F,  G,  H 

1.18 

1.054 

2.054 

3.502 

10 

F,  G,  H, 

I, 

J,  K 

2.18 

1.072 

2.071 

3.531 

11 

ALL  11 

2.28 

1.616 

3.538 

5.203 

12 

B,  0,  I 

2.27 

0.088 

0.006 

0.281 

13 

C,  E,  H 

1.59 

0.161 

0.696 

1.485 

14 

C,  D,  I 

1.96 

0.101 

0.034 

0.296 

15 

I 

I 

I 

.1 


j 

' 

J 

.1 


1 

1 

I 

II 
fl 


-31- 


B  magnetometers.  B  measurements  generally  have  a  larger  V  than  mea- 
x  y 

surements,  which  must  result  from  the  fact  that  the  row  of  magnetometers 

is  displaced  from  the  dipole  in  the  y  direction.  Identification  of  the 

trends  in  the  table  can  be  simplified  in  Fig.  A. 5  by  plotting  V  for  each  type  of 

-*► 

measurement  (B  ,  B  ,  or  B)  as  a  function  of  the  mean  distance  from  the 
x  y 

dipole  to  the  measurement  locations.  In  order  to  separate  an  increase  in 

V  due  to  a  decrease  in  mean  distance  from  one  due  to  an  increase  in  the 

number  of  magnetometers,  the  data  are  sorted  by  number  of  signals.  For 

example,  two  B  measurements,  two  B  measurements,  and  one  B  measurement 
x  y 

each  have  two  signals.  The  graph  shows  several  trends.  For  any  series 

of  measurements,  V  decreases  with  increasing  mean  range;  the  line  through 

points  2v,  3v,  and  4v  has  a  slope  of  -6,  consistent  with  Eq.  4.23  for  a 

single  vector  magnetometer.  Except  at  large  mean  distances,  VD  >  V„  . 

By  Bx 

Whether  this  is  true  in  general  may  depend  on  the  choice  of  measurement 
locations.  The  graph  also  shows  that  high  V  can  be  obtained  at  large  mean 
distances  by  using  a  large  number  of  measurements.  However,  V  is  5.144 

for  five  B  measurements  close  to  the  dipole,  whereas  adding  the  six  more 

*► 

distant  ones  increases  V  only  by  1%.  For  the  3  closest  B  measurements,  V 
equals  4.445,  only  14%  below  the  value  for  5  vector  measurements. 

The  E  and  G  sensitivity  vectors  for  a  B^  measurement  were  shown  to 
be  nearly  perpendicular  in  Fig.  4.4.  A  simple  quadradic  equation  can  be 
solved  to  show  that  the  B^  vectors  will  be  perpendicular  for  two  pairs  of 

measurement  locations:  x  -  ±  0.56y  and  x  ■  ±  3.56  y.  At  the  first,  which 

% 

is  close  enough  to  be  useful,  R  is  diagonal  with  equal  eigenvalues  and 

V  ■  1.44.  This  is  plotted  in  Fig.  4.5  as  point  16y.  However,  the  choice 
of  this  pair  of  locations  over  E  and  G  would  depend  on  the  relative  impor- 

'V 

tance  of  R  being  diagonal  versus  V  being  20%  larger.  The  fact  that  B 


■  ■ 

■  - 


-33- 


measurements  F  and  G  (point  9y)  and  measurements  E  and  G  (point  5y)  have 
almost  the  same  V  while  having  x-direction  separations  of  0.5  and  1.0, 
respectively,  indicates  that  the  apparent  linear  relation  between  log  V 
and  log  r  doesn't  necessarily  hold  for  changes  in  measurement  configura¬ 
tion  in  the  near  field  region. 

In  Fig.  4.5,  the  points  5v,  9v,  and  17v  are  for  two  vector  magneto¬ 
meters,  with  17v  corresponding  to  both  magnetometers  located  at  the  same 
point  (0.0,  1.0),  for  a  V  of  4.0.  This  indicates  that  while  a  large  V  is 
desired,  it  is  possible  to  artificially  enlarge  V  by  repeating  a  measure¬ 
ment  twice.  This  could  be  avoided  by  dividing  V  by  the  number  of  measure¬ 
ments  taken,  which  would  show  B  and  B  at  position  F  to  have  a  V/n  of  1.0 
’  x  y 

and  B  at  positions  E,  F  and  G  to  have  a  V/n  of  0.99. 


4.4  The  Effects  of  Noise 

In  the  preceding  sections,  we  analyzed  several  different  examples  and 
studied  the  variation  in  V  as  a  function  of  distance.  If  the  measurements 
are  affected  by  either  external  or  instrument  noise,  then  the  signal  at 
large  distances  will  be  masked  by  noise,  so  that  such  a  measurement  can  not 
contribute  fully  to  the  modelling  process.  This  effect  can  be  accounted 
for  by  defining  the  total  information  matrix  by 


^TOT 


l  w  R 
i 


(4.24) 


where  wi  is  a  weighting  function  that  is  zero  for  signals  buried  in  noise 
and  1  for  clean  signals.  For  example,  we  can  define  a  suitable  w  by 

Si  ~  nl  ni 

v±  ”  — - -  -  1  -  —  for  nA  <  s1  (4.25) 

w^  “  0  for  n^  >  Sj^ 

where  s^  and  n^  are  the  signal  and  noise  amplitude  for  the  i^ 
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measurement.  Since  the  magnitude  of  the  signal  is  determined  by  both  the 
strength  of  the  source  and  the  sensitivity  vector  for  that  measurement,  we 
can  replace  n^/s^  by  | | / | | ,  where  N^  is  a  noise  sensitivity  vector 
that  contains  the  relative  strength  of  the  noise  source  to  the  signal 
source.  The  use  of  | |  is  reasonable,  since  we  showed  previously  that 
the  magnitude  of  the  sensitivity  vector  indicated  the  maximum  signal  from 
a  dipole  with  unknown  orientation.  From  Eq.  (4.22),  we  find  that,  for  a 
single  vector  magnetometer 


% 

R 


Inb  I 


(4.26) 


It  follows  from  Eq.  (3.8)  that 


lSB  t  -  ~  Jax2  +  y2  |SB  |  -  -J-  Jx2  +  4y2 

x  r  *  y  r  * 


(4.27) 


Suppose  that;  for  a  given  dipolar  source,  the  signal  to  noise  ratio  for  a 
measurement  of  on  the  y  axis  is  unity  for  y  =  2.0.  At  that  point 


l;,B  I 

X 

|sR  I 


so  In, 


0.125 


(4.28) 


-r  -► 

If  we  assume  that  |n  |  ■  |N_  |,  we  can  calculate  how  V  for  the  vector  mea- 

y  “x 

surement  depends  on  distance  along  the  y  axis.  The  results,  plotted  in 
Fig.  4.6,  show  that  at  y  -  1.0,  V  corrected  for  noise  is  10%  smaller  than 
the  original  value.  At  y  -  1.5,  V  is  32%  smaller.  At  y  -  0.5,  the  correc¬ 
tion  is  only  1%. 
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As  a  final  example  of  the  addition  of  noise  to  the  analysis,  we 

stated  earlier  that  V  for  a  single  vector  magnetometer  depended  only  on 

distance  from  the  dipole.  If  the  noise  field  is  not  the  same  for  both 

B  and  B  ,  i.e.  |n^  |  >  |n^  |,  then  there  will  be  an  optimum  magnetometer 
X  y  x  y 

location.  Continuing  the  previous  example,  let  us  assume  that  |N  |  =  0.250 

D 

,  x 
and  |N  |  =  0.125.  On  the  x-axis,  the  total  information  matrix  is  given  by 
B 

y 


4(1- 


x3) 


% 

R  * 


(1-  Ng  X3), 

y 


(4.29) 


while  that  on  the  y-axis  is 


I 

I 

i: 

i 

i 

I 

l 

[ 

I 


(4.30) 


At  the  points  (1,  0)  and  (0,  1),  V  has  the  values  1.75  and  1.68,  respec¬ 
tively,  compared  to  1.81  when  both  noise  sensitivities  were  0.125.  This 
is  consistent  with  Fig.  3.2,  since  the  larger  noise  in  B^  can  be  partially 
offset  by  measuring  Bx  where  it  is  the  largest  —  i.e.  on  the  x-axis. 

If  the  weighted  information  matrix  is  used  in  analyzing  the  row  of 
magnetometers,  the  curves  in  Fig.  4.5  will  be  seen  to  fall  off  more 
quickly  with  mean  distance.  A  detailed  analyf i«;  of  various  magnetometer 
configurations  will  require  realistic  estimates  of  signal  and  noise 
strengths,  to  insure  that  |N^|  is  specified  correctly.  Comparison  of 


'  7"  ‘ 

•  *  **  'if- ■ 


expected  performance  with  that  observed  from  actual  field  measurements 
might  allow  the  use  of  an  empirical  noise  weighting  function,  rather  than 
given  by  Eq.  (4.25).  Thus  we  have  developed  techniques  for  using  sensi¬ 
tivity  vectors  for  analyzing  and  comparing  various  magnetometer  configura¬ 
tions,  both  in  the  ideal  noise-free  case,  and  when  noise  limits  the 
distance  at  which  measurements  can  be  made.  Our  next  step  is  to  relate 
this  work  to  the  reciprocity  theorem  prior  to  beginning  quadrupole 


analysis. 


-38- 


V.  QUADRUPOLES  AND  HIGHER  MOMENTS 

The  preceeding  analysis  used  the  sensitivity  vector  obtained  by  com¬ 
puting  the  gradient  of  the  field  with  respect  to  the  model  parameters. 

The  itl*  measurement  was  shown  to  be  given  by 

'Xj  o’  'Xj 

v±  =  M  (5.1) 

For  a  dipole  model,  vector  notation  lets  us  write  this  as 

P1  -  S£  •  M  (5.2) 

However,  Figures  4.3  and  4.4  showed  that  the  components  of  these  sensiti¬ 
vity  vectors  have  the  same  spatial  dependence  as  the  field  of  a  dipole. 
This  and  Eq.  (5.2)  are  simply  demonstrations  of  the  reciprocity  theorem, 
which  will  be  shown  to  provide  the  basis  for  interpreting  the  sensitivity 
of  a  magnetometer  to  magnetic  quadrupoles  and  higher  order  moments. 


For  a  source  which  is  a  distributed  magnetization,  we  can  identify  the 
magnetization  M(r)  as  a  magnetic  dipole  density  and  integrate  Eq.  (5.3) 
over  the  source  volume  to  obtain 
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Alternatively ,  the  field  can  be  determined  by  noting  that 


B(r)  =  -V<Kr) 


(5.5) 


with  <p  being  described  by  the  traceless  tensor  multipole  expansion  listed 


in  Table  5.1.  The  separation  of  the  source  into  dipole,  quadrupole  and 
higher  moments  is  ideal  for  modelling,  and  we  thus  see  the  need  for  expres¬ 


sions  describing  the  sensitivity  vectors  that  relate  a  field  measurement 
to  each  multipole  moment.  Based  on  our  previous  discussion,  we  could 
compute  the  gradient  of  each  potential  term  in  Table  5.1  to  obtain  the 
corresponding  contribution  to  B,  and  then  compute  the  gradient  of  each  B  with 

respect  to  the  moment  to  obtain  the  sensitivity  vector.  Rather  than 

— ► 

do  this,  we  will  use  the  magnetic  vector  potential  A  to  derive  an  ex¬ 
pression  equivalent  to  Eq.  (5.4)  This  expression  will  help  explain  the 
physical  significance  of  sensitivity  vectors. 

If  we  represent  a  single  axis  magnetometer  by  a  small  pick-up  loop,  then 
the  magnetometer  output  will  be  proportional  to  the  flux  <j>  coupling  the 


magnetometer,  given  by 


/b(t)  • 


(5.6) 


where  da  is  the  normal  to  an  element  of  surface  bounded  by  the  loop. 
-► 

Because  7  •  B  *  0,  we  can  write 


B(r)  -  7  x  A(r) 


(5.7) 


so  that 


*  -  /  7  x  A(r)  •  da 


(5.8) 


Stoke' 8  Theorem  can  be  used  to  convert  this  to  the  line  integral 


*  ■  f  A(r)  •  di 


(5.9) 


where  dl  is  a  length  element  of  the  pick-up  coil  system.  We  now  need  to 


■MfeMpg 
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TABLE  5.1 


I 


The  Traceless  Tensor  Multipole  Expansion  for  a  Distributed  Magnetization 

3  3  3 

Hr)  =  ^  Vi  +  E  E  Qij*ij  +  *  •  * 
i-1  i-1  J “I 

DIPOLE 
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relate  A  to  II,  and  do  this  by  starting  with  Maxwell's  fourth  equation 


which  becomes 


V  x  B  =  p  V  x  M 
o 


V  x  (V  x  A)  =  pQ(V  x  M) 


(5.10) 


(5.11) 


But  the  vector  identity  7  x  (7  x  A)  =  7(7  •  A)  -  V2A  lets  us  write 


V2A  =  V(V  •  A)  -  uq(V  x  M) 


(5.12) 


The  Coulomb  gauge  can  be  chosen  to  set  the  first  term  on  the  right  to 


zero,  leaving 


V2A(r )  =  -\i  (V  x  M(r) ) 
0 


(5.13) 


which  is  Poisson's  equation.  In  cartesian  coordinates,  the  solution  is 


Po  f  V'  x  M ( r ' ' 

'T'j  IJ-M 


(5.14) 


The  vector  identity  VV  x  u  »  V  x  (VU)  -  VV  *  U  can  be  used  to  rewrite  this 


"J  I?-?1!  4"  J  \\f-P\) 


x  M ( r ' ) d 3 r  (5.15) 


The  first  volume  integral  can  be  converted  into  a  surface  integral  that 
vanishes  for  any  bounded  magnetization  distribution.  Therefore 


)  X  7' 


-  r'| 


(5.16) 


, 
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Substituting  this  into  Eq.  (5.9),  we  find 


4rr 


// 


M(r')  x  V 'f  —  W3r' 


r  -  r' 


di 


(5.17) 


The  terms  of  the  scalar  triple  product  can  be  permuted  and  the  order  of 
integration  reversed  to  obtain 


$ 


)  dJr 


3^< 


(5.18) 


Recall  that  a  test  current  1^  flowing  in  the  pick-up  coil  has  a  vector 


potential 


.(?•) 


di 


4ir  7  (-*• 

7  I r  -  r 


(5.19) 


associated  with  a  magnetic  field 


B  (r’)  -  V’x  A  (r‘) 
P  P 


(5.20) 


With  this  identification,  we  may  write  Eq.  (5.18)  as 


t  «=  —  /V’)  *  M(r’)  d3r ’ 


(5.21) 


Thus  we  see  that  for  an  arbitrary  pick-up  loop  configuration,  the  contri¬ 
bution  to  the  magnetometer  output  from  a  source  element  at  can  be 
determined  from  the  magnetic  field  produced  at  r’  by  passing  a  test 
current  1^  through  the  pick-up  coil.  The  reversal  of  the  order  of  inte¬ 
gration  prior  to  Eq.  (5.18)  and  the  identification  of  the  line  integral 
in  Eq.  (5.19)  are  simply  a  statement  of  the  electromagnetic  reciprocity 
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theorem.  Comparison  of  Eqs.  (5.2)  and  (5.21)  results  in  this  identifica¬ 
tion: 

At  a  given  source  location,  the  sensitivity  vector  for  a 
single  axis  magnetometer  is  equivalent  to  the  magnetic 
field  produced  at  that  location  by  passing  a  unit  test 
current  through  the  magnetometer  pick-up  coil. 

Since  the  magnetometers  used  in  Fig.  4.3  and  4.4  were  assumed  to  be  small, 

it  is  correct  that  their  field  B  would  be  dipolar.  Interpretation  of 

P 

Fig.  3.1  is  also  straightforward.  The  sensitivity  vector  for  each  Bx  mag¬ 
netometer  location  corresponds  to  the  magnetic  field  produced  at  the 
origin  by  passing  a  test  current  through  the  magnetometer. 

Equation  (5.21),  with  its  simple  interpretation,  provides  the  key  for 
understanding  the  sensitivity  to  higher  moments:  the  integrand 


B  (r  )  •  M(r  )  can  be  broken  down  into  various  multipole  terms  using  a 
P 

-*■  -*•  -+ 

Taylor's  expansion  of  Bp(r')  and  a  multipole  expansion  of  H(r'),  so  that 

-*■  -+ 

Bp  represents  the  sensitivity  to  the  dipole  component  m,  each  spatial  gra- 
dient  of  B^  represents  the  sensitivity  to  the  corresponding  quadrupole 
moments  and  so  forth.  If  the  magnetometer  is  a  simple  scalar  instrument, 


Bp  can  be  obtained  from  Eq.  (2.3),  where  is  the  unit  dipole  moment 

produced  by  I  .  For  a  B^  magnetometer,  m  =  i.  For  a  gradient  magneto- 
->■ 

meter,  Bp(r)  will  be  characteristic  of  a  magnetic  quadrupole,  but  the 
sensitivity  of  this  instrument  to  a  magnetic  dipole  will  still  be  given  by 


B  •  m. 

P 

Based  on  this  discussion,  we  will  proceed  to  determine  the  sensitivity 
of  a  vector  magnetometer  to  a  dipole-plus-quadrupole  source.  Calculation 
of  the  components  of  Bp  is  straightforward.  There  are  nine  gradients  of 
Bp  that  can  be  written  as  a  gradient  tensor 
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3z 

3z 

3z 

(5.22) 


Just  as  the  quadrupole  tensor  in  Table  5.1  is  traceless  and  symmetric,  the 
— ► 

condition  that  V  *  B  is  zero  in  a  current-free  region  implies  that  the 

-y  -*• 

gradient  tensor  for  Bp  is  symmetric,  so  that  we  can  write  the  B^  gradient 

tensor  as 


3B 

_ X 
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3B 

3x 
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(5.23) 


The  condition  that  V  •  B^  =  0  implies  that  this  tensor  is  also  traceless. 

We  can  therefore  omit  both  3B  /3z  and  Q  from  our  discussion,  since  they 

z  z  z 

can  be  derived  from  the  other  two  terms  of  each  trace. 

-^L’he  following  pairs  of  quantities  will  be  required  for  the  sensitivity 
analysis 


1: 

B  , 

m 

X 

X 

OIPOLE  ^ 

2: 

B  , 
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y 

(5.24) 
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(5.24  cont'd) 


4:  — -  .  Q 

9x  ’  wxx 


5:  aT  •  Qyy 


QUADRUPOLE 


6:  9x  ’  Qxy 


71  JT  ’  Qxz 


8:  —  ,  Q 
9y  y* 


The  terms  on  the  right  correspond  to  the  eight  axes  in  model  space,  while 
the  quantities  on  the  left  correspond  to  the  eight  components  of  each 
sensitivity  vector.  If  each  set  of  eight  quantities  is  thought  of  as  a 
(1  x  8)  matrix  or  an  8-vector,  then  the  magnetometer  output  corresponds  to  the 
inner  (scalar)  product  of  these  two  matrices  (vectors).  In  computing  the 
gradients  of  using  Eq.  (2.3),  it  is  important  to  note  that  r'  is  the 
location  of  the  pick-up  coil  and  r  is  the  location  of  the  model,  and  the  gra¬ 
dients  are  computed  with  respect  to  the  source  location.  The  sign  of  I 
(or  nip)  can  be  chosen  to  insure,  for  example,  that  as  shown  in  Fig. 

5.1  will  produce  a  positive  Bx  on  the  positive  x-axis. 

The  eight  sensitivity  vector  components  for  each  of  the  vector  mag¬ 
netometer  axes  are  listed  in  Table  5.2.  These  results,  which  were 
obtained  by  computing  the  appropriate  gradients  of  B^ ,  could  also  be  ob¬ 
tained  by  computing  the  appropriate  gradients  of  the  potentials  in  Table 

5.1.  Comparison  of  the  two  methods  shows  that  a  factor  of  p  /4it  must 

o 

multiply  each  term  in  Table  5.2  to  provide  the  proper  units,  and  an 
additional  factor  of  1/2  must  multiply  the  quadrupole  terms  to  provide  a 


- 

.  ....  V  . 


*  ■ 


TABLE  5.2 


The  sensitivity  vector  components  for  each  term  of  a  dipole-quadrupole 

model  and  each  axis  of  a  vector  magnetometer.  Each  term  should  be  multi 
U 

plied  by  - —  ;  the  quadrupole  terms  should  be  multiplied  by  a  factor  of  *s 
4tT 


TABLE  5.3 

The  sensitivity  vector  components  for  each  axis  of  a  vector  magnetometer 
located  at  the  points  (1,  0,  0),  (0,  1,  0)  and  (0,  0,  1)  obtained  from 
the  equations  in  Table  5.2. 


(1,  0,  0) 


(0,  1,  0) 


0 

0 

2 

0 

0 

-1 

-6 

0 
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0 

0 

-6 

(0,  0,  1) 
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correct  Taylor's  expansion  of  Bp. 

Table  5.3  lists  the  sensitivity  vector  components  for  vector  magneto¬ 
meters  a  unit  distance  from  the  source  along  the  positive  x,  y,  and  z 
axes.  It  is  important  to  note  from  the  table  that  Q  does  not  affect 
B(l,  0,  0),  while  Q^z  and  Q^z  do  not  affect  B(0,  1,  0)  and  B(0,  0,  1), 
respectively.  This  is  in  contrast  to  the  single  dipole  case  discussed 


previously,  in  which  the  invertability  of  the  field  equations  was  depen¬ 


dent  only  upon  distance  from  the  source.  Table  5.3  can  also  be  used  to 
clarify  the  orientation  of  the  field  from  each  quadrupole  in  Fig.  5.1. 


Given  the  sensitivity  vector  equations  in  Table  5.2,  the  previous 
numerical  examples  can  be  readily  extended  to  include  the  quadrupole  terms. 
However,  the  information  matrix  will  then  be  (8x8),  which  can  be  manipu¬ 
lated  most  efficiently  with  a  digital  computer  and  thus  will  not  be 
addressed  further  in  this  report. 
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VI .  NONLINEAR  PROBLEMS 

In  all  of  the  models  discussed  so  far,  the  field  equations  were 
linear  with  respect  to  the  source  parameters.  For  this  reason,  the  com¬ 
putation  of  the  sensitivity  vectors  was  simple,  and  the  results  were 
consistent  with  the  reciprocity  theorem.  If  we  had  chosen  the  position 
of  the  dipole  as  a  model  parameter,  the  three  components  of  r'  in  Eq.  (2.3) 
would  represent  three  non-linear  terms  in  the  field  equation.  We  will  now 
proceed  to  extend  our  sensitivity  vector  techniques  to  include  this  case. 


6. 1  Taylor's  Series  Linearization 

In  section  II,  we  showed  that  for  linear  systems  the  fields  and  model 
are  related  by  an  equation  of  the  form 

a.  >v  ^ 

F  =  T  M  (6.1) 

y  ^  a. 

where  F  is  an  n  x  1  matrix,  M  is  an  m  x  1  matrix,  and  T  is  m  x  n.  In  the 

non-linear  case,  we  can  write 

'V  'Xj  'V 

F  =  f(M)  (6.2) 

%  'V, 

where  f(M)  represents  a  set  of  n  non-linear  functions  of  the  m  model  para- 

meters.  For  example,  each  equation  in  f  could  be  a  different  polynomial. 

If  we  choose  a  point  in  model  space  M  and  determine  the  corresponding 

^o 

point  in  field  space  F 

%o  ^ 

F°  =  f(M°)  (6.3) 


If  we  restrict  ourselves  to  the  point  M  ,  we  can  write  this  as 

F°  “  T°  M°  (6.4) 

We  can  then  approximate  f  in  a  region  about  this  point  using  a  Taylor's 
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series  expansion  of  each  of  the  n  equations  in  f.  In  matrix  notation. 


these  can  be  written  as 


f\j  f\,  % 

F  -  f (M) 


A.)  ♦ 


(„  .  + 

j  *]  1  i 


(6.5) 


The  summation  containing  the  partial  derivative  can  be  represented  by  a 
matrix  multiplication,  so  that  the  equation  reduces  to 


>\,  %  f\j  *\j  *\j  *\j 

F  *  F°  +  J(M°)(M  -  M  )  +  .  . 


(6.6) 


where  J(M°)  is  the  Jacobian  matrix  J(M)  evaluated  at  M  ,  with 


•  '  8M 


(6.7) 


The  right  hand  side  of  Eq.  (6.5)  is  a  linear  vector  function  of  M  which 
best  approximates  the  non-linear  function  f(M)  at  the  point  M  .  This 

'Xj  % 

linear  function  describes  the  tangent  hyperplane  to  the  surface  F  =  f(M) 


If  we  know  F  from  a  series  of  measurements  and  wish  to  determine  the 

corresponding  M,  it  is  necessary  to  solve  Eq.  (6.2).  One  way  to  do  this 

^  ^o  V> 

is  make  a  trial  guess  for  M,  say  M  ,  and  use  Eq.  (6.3)  to  compute  F  . 

The  error  in  the  fields  can  be  written  as 


AF  -  F  -  F° 


(6.8) 
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Frora  Eq.  (6.6),  the  error  in  the  model  is  given  by 


'v  %  ^ 

AM  =  M  -  M  =  J(M  )  AF 


(6.9) 


The  choice  of  the  point  for  computing  the  Taylor's  series  can  be  changed 
to  MJ  =  M  +  AM  and  the  process  repeated  iteratively  to  obtain  the  desired 
accuracy.  Of  course,  problems  will  occur  if  the  Jacobian  matrix  is  non- 
invertable  or  if  the  choice  of  M  leads  to  a  phantom  solution. 


6.2  Sensitivity  Vectors  for  Non-Linear  Systems 

If  we  compare  Eq.  (6.5)  with  Eq.  (2.7)  which  defined  sensitivity  vec¬ 
tors,  we  see  that 


m  9f . 

£"7 


(MJ  '  V  + 


(6.10) 


ul  r\j 

ri°  +  L  <si>j  -  “3°> 


(6.11) 


which  by  Eqs.  6.8  and  6.9  becomes 


AFi  =  Sj[  AM 


(6.12) 


Using  Eq.  (6.4)  instead  of  F°,  we  could  also  have  written  Eq.  (6.11)  as 


ri  ■  £  V  »i°  ♦  V  «S  -  «°> 


(6.13) 


It  is  apparent  that  the  sensitivity  vectors  no  longer  correspond  to  the 

'V. 

rows  of  T,  but  instead  to  the  rows  of  the  Jacobian  matrix. 

The  derivations  of  Section  IV  can  be  repeated  to  relate  this  to  the 


a. 


■ 

' 
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information  matrix.  Multiplying  Eq.  6.12  by  and  summing  oyer  all  mea¬ 
surements,  we  obtain 


L 


'v*  *  T  ~ 

S±  S±  AM 


(6.14) 


i-1 


i-1 


This  can  be  solved  for  AM 


AM 


n  V1  A 

Y  S  s  T  )  Ys±  AF 
i-1  1  1  /  i-1 


(6.15) 


Using  Eq.  (6.11),  we  can  also  write 

n 


'V*  'V>o 
M  =  M 


R"1 


(6.16) 


i-1 

This  equation,  which  uses  sensitivity  vectors  and  the  information  matrix, 
represents  the  solution  of  the  linearized  version  of  our  original  non¬ 
linear  field  equation. 

Thus  we  see  that  even  with  non-linear  systems,  sensitivity  vectors 
have  a  straightforward  mathematical  interpretation.  As  in  the  linear  case, 
their  utility  lies  in  their  simple  physical  interpretation  in  terms  of 
measurement  sensitivity  and  in  their  suitability  for  graphical  presenta¬ 
tion.  It  is  important  to  recognize  that  in  practice,  sensitivity  vectors 
can  be  determined  either  analytically  or  numerically.  If  we  have  a 
measurement  that  is  a  linear  function  of  some  variables  and  a  non-linear 
function  of  several  others,  we  can  use  the  sensitivity  vector  approach  to 
determine  how  the  measurement  is  affected  by  each  variable,  regardless  of 
whether  the  variable  has  a  linear  or  non-linear  effect  on  the  measurement. 
However,  we  have  just  shown  that  the  interpretation  does  depend  somewhat 
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on  whether  the  dependence  is  linear  or  non-linear.  If  the  i*"^  measurement 
is  described  by  F^,  it  follows  from  our  original  definitions  that  for  the 
linear  case ,  the  sensitivity  vector  for  this  measurement  is 


Si  "  Fi 


(6.17) 


and  that 


Fi'  Si 


-*■ 

M 


(6.18) 


Because  of  the  linear  relationship,  the  change  in  the  measurement  result¬ 
ing  from  a  change  in  the  model  Is  simply 


AFi  =  Si  •  AM 


(6.19) 


In  the  non-linear  case,  the  sensitivity  vectors  is  still  given  by  Eq . 
(6.17),  but  we  have  shown  that  Eq.  (6.18)  does  not  apply,  and  that  (6.19) 
does  apply  as  long  as  higher  order  terms  can  be  neglected.  Thus  we  can 
write 

non-linear  case,^ 

-► 

small  AM 


AFi  *  Si 


AM 


(6.20) 


This  leads  us  to  the  important  conclusion 

If  sensitivity  vectors  are  interpreted  as  a  measure 
of  how  a  measurement  is  affected  by  a  small  change 
in  the  model,  then  linear  and  non-linear  problems 
can  be  treated  identically. 


6.3  The  Sensitivity  of  a  Vector  Magnetometer  to  Dipole  Position 

As  a  non-linear  example  of  sensitivity  vectors,  we  will  extend  our 
previous  example  and  determine  the  sensitivity  of  a  vector  magnetometer 
to  small  changes  in  the  location  of  the  dipole.  The  sensitivity  of  Bx  to 
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changes  in  dipole  location  can  be  determined  by  computing  the  gradient  of 
Bx(r),  given  by  Eq.  (2.4),  with  respect  to  r  * .  The  resulting  sensitivity 
vector  is 


5<V 


SB  (r) 

X 

3x* 


- 1 
i 


3B  (r)  . 
”  r-  j' 


3y 


(6.21) 


This  can  be  combined  with  the  sensitivity  vector  for  the  dipole  components, 
to  obtain 


S(B  ) 

X 


3V;> 

3m 


3B  (r)  „ 

+  -* -  j 

3m  J 


3B  (r)  „  3B  (r)  „ 

+  -4i-  i'  +  4,  j’ 


3x' 


3y' 


St  +  S_j  +  S-i'  +  S , j ' 


(6.22) 


The  analytic  expressions  for  these  sensitivity  vector  components,  as  well 
as  the  corresponding  ones  for  B^,  are  listed  in  Table  6.1.  It  is  impor- 
tant  to  note  from  Eq.  (2.3)  that  the  gradient  of  B  with  respect  to  r  is 
equal,  with  opposite  sign,  to  the  gradient  with  respect  to  r'.  Thus 
and  also  describe,  if  multipled  by  -1,  the  change  in  the  measurement 
produced  by  a  small  change  in  the  magnetometer  location. 

The  equations  for  and  indicate  that  the  position  sensitivity 
at  a  particular  location  depends  upon  the  dipole  components.  For  this 
reason,  we  must  choose  in  for  example  of  position  sensitivity.  As  an  ex¬ 
tension  of  Table  3.3  and  Figure  3.2,  Table  6.2  lists  B  ,  B  ,  the  corres- 

x  y 

ponding  values  of  and  S^,  and  the  quantity  for  2  dipole  orientations 
where 

V34  "  lS3CBx>S4CV  '  WWl  (6.23) 

Comparison  of  the  m  *  (1,  0)  data  for  locations  A  and  E  shows  that  a  vector 
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TABLE  6.1 


I 

! 

[ 

( 


The  four  components  of  the  sensitivity  vectors  of  Bx  and  By  magnetometers 

that  describe  how  B  and  B  are  affected  by  small  changes  in  the  components 
x  y 

and  position  of  a  dipole  located  at  the  origin. 


S1(BX) 


I 

I 


I 

( 


9B 

W  *  ir  -  -  (3x*> 

y  r 

3B  1 

S3(V  =  =  "V  [mx(6x3 " 9xy2)  +  my(12x2y  -  3y3)J 

s4(bx)  ”  jp-  ■  ~  [mx(12x2y  -  3y3>  +  my(-3x  +  12xy2>J 

9B  . 

si<V  ‘  S?  '  -T  <*"> 

7  x  r 

9B  . 

W  '  ‘  "7  -  r2> 

y  r 

3B  *1 

s3(By)  "  9^  “  “7  [mx<12x2y  -  3y3>  +  my(_3x3  +  12xy2)J 


position  sensitivity  vectors  So  and  S,  of  B  and  B  measurements  for  two  dipole  orientations 
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field  measurement  is  twice  as  sensitive  (V.^  =  18)  to  changes  in  dipole 
position  for  measurements  on  the  dipole  axis  than  for  equidistance  mea¬ 
surements  perpendicular  to  the  dipole  axis  (V,^  =  ®)»  The  position 
sensitivity  is  even  lower  at  point  D  (V_^  =  6.0).  The  relationship  of 
the  position  sensitivity  vectors  for  Bx  and  measurements  is  shown 
graphically  in  Fig.  6.1. 

While  and  follow  directly  from  the  slope  of  vs  x,  Bx  vs  y. 

By  vs  x  and  B^  vs  y  signatures  passing  through  each  measurement  point,  the 
orientation  of  the  position  sensitivity  vectors  follows  more  directly  from 
another  type  of  plot.  From  Eq.  (6.20),  the  change  in  field  associated 
with  a  change  in  position  will  be 


ABi  =  Si  •  AM^  =  S3Ax'  +  S^Ay' 


V'  Bi(r)  •  Ar' 


(6.24) 


VB±  •  Ar 


Thus  any  displacement  of  the  magnetometer  parallel  to  a  line  of  constant 
B^  will  produce  no  change  in  field.  A  corresponding  displacement  in  the 
source  will  also  produce  no  change.  By  identifying  the  sensitivity  vec¬ 


tor  with  the  spatial  gradient  of  B^,  we  see  immediately  that  must 
be  perpendicular  to  lines  of  constant  B^.  Given  that 


p  m(2x2  -  y2) 

Bx  -  4ir  5 

r 


u  3mxy 

B  *  - - 

y  4tt  r5 


(6.25) 


it  follows  from  the  parametric  substitution  y  =  ax  that  the  lines  of  con¬ 
stant  Bx  must  satisfy  the  equations 


1  m(2  -  a2) 

B  4ir  c  /  o 

x  ,,  9.5/ 2 

a  +  a2)  , 


gjnirapvsi - - 

■ 


y  •  ax 


(6.26) 


Fig.  6.1  The  position  sensitivity  vectors  for  measurements  of  B  and  B 


at  unit  distance  from  a  magnetic  dipole  located  near  the  origin  and  par 


allel  to  the  x-axis.  The  data  are  obtained  from  Table  6.2 
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while  those  for  B^  must  satisfy 


1  p  m(3a) 
_ o  _ 

B  4tt  (1  +  a2) 

y 


y  =  ax 


(6.27) 


The  lines  where  the  magnitude  of  B  is  constant  satisfy 


B  4tt 


4  +  a' 


(1  +  a2)1 


,  y  =  ax 


(6.28) 


One  set  of  iso-B^  lines  is  plotted  in  Figure  6.2,  along  with  the  corres¬ 
ponding  B^  sensitivity  vectors  for  the  points  in  Fig.  6.1.  Clearly,  the 
position  sensitivity  vector  corresponds  to  the  gradient  vector  at  each 
point.  Thus  a  complete  set  of  iso-B^  curves,  in  the  form  of  a  contour 
map,  would  be  an  efficient  means  of  visually  identifying  the  magnetometer 
location  with  either  high  or  low  position  sensitivity.  Furthermore,  since 
iso-Bx  lines  need  not  be  drawn  for  fields  lower  than  the  ambient  noise, 
this  type  of  plot  will  also  identify  the  zones  where  B^  can  be  measured 
with  adequate  signal  to  noise  ratio. 

The  iso-By  curves  are  shown  in  Fig.  6.3,  with  the  obvious  difference 

that  B  from  an  x  dipole  is  zero  on  the  coordinate  axis  while  B  is  a 
y  x 

~¥- 

maximum  there.  Figure  6.4  shows  the  iso-|B|  curve,  which  can  be  seen  to 
bound  the  curves  for  B^  and  B^.  This  figure  shows  that  an  instrument  that 
measures  |b|  directly  would  have  a  much  simpler  position  response  than  one 
measuring  only  a  single  component.  It  is  also  important  to  note  that  by 
the  reciprocity  theorem  each  of  these  iso-field  curves  can  also  be  thought 
of  as  an  iso-sensitivity  curve  for  a  magnetometer.  This  concludes  the  ana¬ 
lysis  of  non-linear  sensitivity  vectors.  Because  of  their  straightforward 
interpretation,  the  iso-field  plots  and  sensitivity  vectors  may  be  especially 
useful  for  studying  the  position  sensitivity  of  quadrupole  models. 


hK ; 


VII.  CONCLUSIONS 


In  this  report,  the  mathematical  basis  of  sensitivity  vectors  has 
been  developed,  with  specific  application  to  the  measurement  and  modelling 
of  magnetic  fields.  In  the  course  of  analyzing  several  examples,  the 
following  questions  were  answered: 

1.  Given  two  magnetic  field  measurements,  what  criteria  must 
be  satisfied  to  allow  determination  of  the  two  dipole  com¬ 
ponents  consistent  with  that  field? 

2.  For  a  measurement  of  the  vector  magnetic  field  at  a  fixed 
distance  from  a  dipole  source,  what  is  the  optimum  magne¬ 
tometer  position? 

3.  How  do  various  configurations  of  single  axis  magnetometers  compare 
in  terms  of  their  ability  to  determine  model  parameters, 

and  how  can  this  be  assessed  quantitatively? 

4.  How  can  the  presence  of  noise  be  included  in  such  an 
analysis? 

5.  What  is  the  relationship  of  the  sensitivity  vector  concept 
to  the  reciprocity  theorem  of  electromagnetic  fields? 

6.  How  can  sensitivity  vectors  be  used  to  study  multipole 
models? 

7.  What  is  the  interpretation  of  a  sensitivity  vector  for  a 
non-linear  model? 

8.  What  is  the  position  dependence  of  a  vector  magnetometer, 
and  how  does  it  vary  with  position? 

9.  Where  can  a  magnetometer  be  placed  relative  to  the  model 
to  insure  adequate  signal  to  noise  ratio? 

Based  on  the  ability  of  the  sensitivity  vector  approach  to  answer  these 
questions  quantitatively,  it  appears  that  this  type  of  analysis  may  be 
valuable  for  optimizing  magnetometer  array  configurations.  It  is  re¬ 
assuring  that  the  results  obtained  in  this  report  are  consistent  with  both 
a  modeller's  intuition  and  more  abstract  mathematical  analysis,  particular¬ 
ly  in  that  this  method  can  be  readily  extended  to  more  complicated  systems 


where  intuition  fails. 


Two  specific  recommendations  follow  from  this: 

1)  Interactive  computer  code  should  be  developed  to  allow 
accurate  and  rapid  analysis  of  magnetometer  sensitivity 
and  comparision  of  magnetometer  configurations. 

2)  The  sensitivity  vector  analysis  should  be  extended  to 
include  the  magnetic  field  from  electric  current  dis¬ 
tributions,  and  possibly  the  electric  field  from  these 
currents. 

In  conclusion,  sensitivity  vectors  are  a  potentially  valuable  tech 
nique  that  can  be  readily  applied  to  the  measuring  and  modelling  of  mag 
netic  fields.  Its  greatest  benefit  over  other  techniques  is  that  it 
allows  simple,  graphical  display  of  magnetometer  sensitivity. 
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VIII.  FIGURE  CAPTIONS 


Fig.  3.1  The  sensitivity  vectors  for  measurements  of  B^ 

at  unit  distance  from  a  fixed  magnetic  dipole. 
The  coordinates  of  the  points  are  listed  in 
Table  3.1 . 


Fig.  3.2  The  sensitivity  vectors  for  measurements  of 

and  By  at  unit  distance  from  a  fixed  magnetic 

dipole.  The  coordinates  of  the  points  are 
listed  in  Table  3.2 . 


Four  locations  where  the  Bx  sensitivity  vector 

is  perpendicular  to  that  at  one  particular  lo¬ 
cation  and  has  the  same  magnitude . 


A  row  of  B  or  B  magnetometers  near  a 
x  y 


magnetic  dipole. 


The  quantities  S  and  S  plotted  as  a  function 
x  y 

of  x  for  the  B  magnetometers  in  Fig.  4.2  and 
a  vector  plot  of  S  for  the  same  data . 


The  quantities  and  plotted  as  a  function 

of  x  for  the  B  magnetometers  in  Fig.  4.2  and 
7  -f 

a  vector  plot  of  S  for  the  same  data . 


The  quantity  V  plotted  as  a  function  of  mean 
distance  r  from  the  data  in  Table  4.4.  The 
lines  connect  sets  of  measurements  using  the 
same  type  and  number  of  magnetometers.  .  .  . 


The  quantity  V  plotted  versus  y  for  a  vector 
magnetometer  on  the  y  axis  computed  with  and 

without  the  noise  correction  . 


Fig.  5.1  The  source-sink  pictures  for  the  dipole  and 
quadrupole  in  the  traceless  tensor  represen¬ 
tation . 

Fig,  6.1  The  position  sensitivity  vectors  for  measure¬ 
ments  of  B  and  B  at  unit  distance  from  a 
x  y 

magnetic  dipole  located  near  the  origin  and 
parallel  to  the  x-axls.  The  data  are  obtained 
from  Table  6.2 . 


feet*  «■ 
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Figure  Captions  (cont'd) 

Fig.  6.2  The  lines  where  |B  |  produced  by  mx  is  a 

constant .  61 

Fig.  6.3  The  lines  where  |B  |  produced  by  m ^  is  a 

constant . ^ .  62 

Fig.  6.4  The  line  where  |b|  from  is  a  constant .  63 


IX.  LIST  OF  TABLES 


Table 


Sensitivity  vectors  at  several  points  equi¬ 
distant  from  a  dipole.  The  data  are  plotted 
in  Fig.  3.1 . 


The  quantity  V  * 


lSx  SY 
X1  2 


yl  X2 


computed 


for  pairs  of  Bx  measurement  points  in  the 


first  and  second  quadrants. 


The  coordinates  and  sensitivity  vector 

components  for  B  and  B  measurements  in 
x  y 

the  first  quadrant.  The  data  are  plotted  in 
Fig.  3.2 . 

The  information  matrices  for  each  of  the  mea¬ 
surement  points  and  pairs  of  points  in  Table  3.2. 

Sensitivity  vectors  and  information  matrices  for 
the  line  of  Bx  magnetometers  in  Fig.  4.2.  .  .  . 


Sensitivity  vectors  and  information  matrices  for 
the  line  of  B^  magnetometers  in  Fig.  4.2.  .  .  . 
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X.  LIST  OF  SYMBOLS 


Symbols 


A 

-* 

B 

-► 

B 

P 


Magnetic  vector  potential 
Magnetic  Induction 

Magnetic  induction  from  I  In  pickup  coil 
Electric  field 

Scalar  field  measurement  matrix 

Non-linear  functions  of  M 
Test  current 


,  j,  k» 


,  M 
,  m 


i 

i 


x  1 
x  1 


Identity  matrix 
Cartesian  unit  vectors 
Jacobian  matrix 
Model  matrix 

Magnetic  dipole 

Dimension  of  M 
Dimension  of  F 
Noise  sensitivity  vector 

Quadrupole  Tensor 

Information  matrix 

Field  and  Source  point  distance  vectors 
Sensitivity  vector 

Transfer  matrix 


V 


Model  space  volume 
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List  of  Symbols  (cont'd) 

Symbols 

Noise  weighting  function 


a 


o 


V 


V' 


Iso-field  plot  parameter 
Kronecker  Delta 

Magnetic  flux 
Multipole  unit  potential 


Magnetic 


permeability  of  free  space 


Gradient  with  respect  to  r 
Gradient  with  respect  to  r' 
Gradient  with  respect  to  M 
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